DESY 04-195, CPT-2004/P.107, WUB 04-17 hep-lat/0411022 _v2 

Nov 2004 - Jan 2005 



Scaling tests with dynamical overlap 
and rooted staggered fermions 

Stephan Diirr a 5 and Christian Hoelbling crf 

a DESY Zeuthen, Platanenallee 6, D-15738 Zeuthen, Germany 
6 Institut fiir theoretische Physik, Universitat Bern, CH-3012 Bern, Switzerland 
c Centre de Physique Theorique*, Case 907, CNRS Luminy, F-13288 Marseille Cedex 9, France 
d Bergische Universitat Wuppertal, Gaussstr. 20, D-42119 Wuppertal, Germany 

(N 

^ ■ Abstract 

We present a scaling analysis in the 1-flavor Schwinger model with the full overlap and 
the rooted staggered determinant. In the latter case the chiral and continuum limit of the 
scalar condensate do not commute, while for overlap fermions they do. For the topological 
CN ■ susceptibility a universal continuum limit is suggested, as is for the partition function and 

the Leutwyler-Smilga sum rule. In the heavy-quark force no difference is visible even at 
finite coupling. Finally, a direct comparison between the complete overlap and the rooted 
staggered determinant yields evidence that their ratio is constant up to 0(a 2 ) effects. 
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1 Introduction 



An issue which is of both phenomenological and conceptual relevance is whether it is a valid 
approach to use the staggered action to study QCD with iVy = 2oriVy = 2 + l dynamical quarks. 
This is because the staggered Dirac operator D st leads to 4 degenerate flavors in the continuum, 
and the simulations are thus performed by taking the square root of the staggered determinant 
(plus a quartic root for the dynamical strange quark). While the results of such studies look 
very promising from the phenomenological viewpoint PjJ |2j, a field-theoretic justification of 
the rooting procedure might be hard to find. The origin of the problem is that no one has 
constructed, so far, a local operator which, when raised to the fourth power, reproduces D st . 

It has been shown that the most naive choice results in a non-local operator [Sill], but of 
course a more elaborate construction might settle the issue. 

The evidence in favor of the rooting procedure that comes from NLO staggered chiral 
perturbation theory [3] needs to be backed by numerical checks of the associated predictions 
(some are established [2]), and the analytical thoughts in 6 J involve a number of simplifications. 
It has been shown - both in 2D and in 4D - that staggered eigenvalues on smooth enough 
backgrounds form near-degenerate pairs/quadruples which (apart from a rescaling factor) mimic 
the (non-degenerate) overlap eigenvalues on the same configuration IE] • And the consequence 
that rooted staggered fermions satisfy an approximate index theorem and are in the right 
random-matrix universality class has been explicitly verified 0EIJ. However, all these pieces 
are inconclusive, as they do not say how the rooting procedure should be matched in the valence 
sector. 

In the absence of a strict analytic argument, only a careful scaling study has some conceptual 
power, albeit an asymmetric one. If a continuum limit is found with rooted staggered fermions 
which agrees with another approach which is considered conceptually proof, nothing firm can 
be said (albeit a failure of the staggered framework seems less likely then). On the other hand, 
if the continuum limits disagree, the staggered answer would be in conflict with universality. In 
this note we attempt such a scaling study. The safe approach against which we shall compare 
is the overlap formulation [llj, which, however, is far more demanding in terms of CPU time. 
The theory in which we will work is the massive Nf = l Schwinger model, i.e. just QED in 2D. 

We define the massless overlap operator as [11 J 

D w 

D m = p (l + 75 sign( 75 D w p )) = p (l + ' 
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with D^f p the Wilson operator at negative mass — p, and construct the massive overlap via 
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DZ={l--)D^ + m. (2) 
2p 



The massless operator ([I]) satisfies the Ginsparg- Wilson relation 



D% + l5 D = 0, 75 = 75 (1 - -D) (3) 

P 



which substitutes the continuum chiral symmetry by the lattice chiral symmetry group 

5ip = j 5 ip, 5tp = V>7 5 (4) 

which, in turn, excludes additive mass renormalization and prevents operators in different 
chiral multiplets from mixing. On the lattices considered below, the full Dirac matrix may 
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be kept in memory, and one can use standard linear algebra routines to perform a singular 
value decomposition of the shifted Wilson Dirac operator D_ =U SV\ where U, V are unitary 
matrices and S > is diagonal. The massless overlap operator is then simply given by 

D ov =p(l + UV^) (5) 

and it is straightforward to plug in a mass, see (J2J), and call further library routines to determine 
the complete eigenvalue spectrum. In practice, one may prefer to determine the eigensystem of 
£)W £)W or Q f /ygj^w bvrt this does not bring any change in principle. We use p— 1, which we 
checked, following Ref . ^1] , is an almost optimal choice with respect to locality for /3 > 4. 

The massless staggered operator reads 

^ St = nE^^lfeW^+w - Ul(x-ft)8 x ^y) (6) 

with ^(x) = (— l)^"<M a: ' / ) and the massive operator is simply D^ t = D st + m. What remains of 
the continuum SU(4)a chiral symmetry group in the case of (jHJ) is the abelian 

X (x)^exp(iM-l)£^ X (x) > x(x)^x(x)exp(-i0 A (-l)£,^), (7) 

which, however, still protects the fermions against additive mass renormalization. The parallel 
transporter Up(x) in @ may be replaced by a weighted sum of gauge- covariant paths from x 
to x+fi [15J. As a result, the most unphysical effects in the staggered formulation, the "taste- 
changing" interactions due to highly virtual gluon exchanges [HI], can be considerably reduced. 
The reason is a separation of the relevant low-energy modes from the regularization dependent 
(and wildly fluctuating) high-energy modes and this is why we will speak of "UV- improved" 
or "UV-filtered" staggered quarks. For comparison the same modification will be considered in 
the overlap case [13 0E], too, but there the effect will be much smaller. 

We will be interested in the Schwinger model (SM) with one active flavor, since only there 
the rooting issue exists (in 2D the staggered formulation generates 2 flavors in the continuum), 
but for completeness let us mention the relationship to QCD both for Nf = l and Nf>2. 

In the zero temperature SM with Nf — 1 the scalar condensate Xsca = ('0V') a ^ zero quark 
mass follows from the global axial anomaly and is given by [IB] 

X - (m = 0) =^M = 0.1599... (8) 

e 2iv i ' z 

where 7 = 0.5772... is the Euler constant. Finite temperature effects will reduce Xsca HE], but 
no temperature will be large enough to really make it vanish. In other words, for Nf = l there is 
no chiral phase transition, and thus the situation in the SM is analogous to QCD with Nf = l. 

With Nf>2 and a small mass term the zero temperature SM shows a vague similarity to 
QCD slightly above the phase transition: The Polyakov loop does not vanish, and the chiral 
condensate is almost zero; the system "tries" to break the axial flavor symmetry spontaneously, 
the spectrum shows a gap between the "Schwinger particle" with mass eJ 'Nf jix + 0(m) (the 
analog of the rf) and the N f 2 -l light "Quasi-Goldstones" with mass M n ^ m N f/( N f+ 1 ) gJJ. The 
latter get sterile in the chiral limit (as required by Coleman's theorem [21]), but the important 
point is that, as long as one stays away from the chiral limit, these "pions" dominate the 
correlators between external P- and A^-sources at sufficiently long distance. 
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A peculiarity of any 2D gauge theory is that the fundamental coupling has the dimension of 
a mass. Furthermore, the SM is super-renormalizable and the coupling e does not run, hence 
the theory is not asymptotically free. We use these features in taking the liberty to set the 
scale through e. In the compact lattice formulation, the fundamental charge is related to f3 via 

P = T^5 ' ( 9 ) 

(below we will set a = l) and this means that the physical scale in the full theory depends only 
on j3 and not on the fermion mass m. Of course, other choices would have been possible. One 
might set the scale via the measured slope in the short-distance potential (which is finite in 
2D, since lim r ^ V(r)/r = const in the continuum with a constant independent of Nf). In fact, 
our choice is just this short-distance option, up to cut-off effects. 



2 Eigenvalue spectra 

The plan is to produce a quenched ensemble, to compute all eigenvalues of the massless D ov , D st 
on each configuration and to introduce the dynamical fermions through reweighting [22] • We 
prepared a variety of lattices with fixed physical volume, both at zero temperature (i.e. with 
a time-extent T = L large compared to all correlation lengths) and at a temperature given via 
T = L/4, see Tables El for details. In 2D standard APE-smearing already involves the full 
hypercube, and we give the staple and the original link weight 0.5 each. Technically, we thus 
apply a smearing step and evaluate the fermion matrix on the resulting background, but we 
consider it a modification of the fermion action. What is peculiar to the choice of setting the 
scale through \fft is that our lattices are still (exactly) matched after reweighting to Nf = 1,2. 

Fig-E]shows the low-energy spectrum of D st and the A of D ov [cf. eqn. (jllSjl ] on four selected 
configurations at (3 = 7.2. We define their topological charge as the overlap index |2"3"1 IT3*] 

q(U) = md(U) = ^tr( 75J D ov ) • (10) 

The first three are typical for topological charge q = 0, 1,4, respectively, while the last panel 
shows one of the rare (at (3 = 7.2) cases where the overlap charge (jlO)) depends on the filtering 
level or, likewise, on the parameter p. Without smearing there is a vague similarity between the 
staggered and the overlap spectrum on the q = configuration, but not on those with higher 
charge. This holds for fairly smooth gauge fields; the plaquette is 0.927722(54) at (3 = 7.2. 
However, after a single smearing step, the situation improves dramatically. The staggered 



geometry 


8x8 


12x12 


16x16 


20x20 


24x24 


28x28 


P 


0.8 


1.8 


3.2 


5.0 


7.2 


9.8 


statistics 


10,000 


10,000 


10,000 


10,000 


10,000 


10,000 



Table 1: Survey of matched zero temperature lattices, with statistics for both types of fermions. 



geom. 


8x2 


16x4 


24x6 


32x8 


40x10 


48x12 


56x14 


64x16 


72x18 


80x20 


P 


0.8 


3.2 


7.2 


12.8 


20.0 


28.8 


39.2 


51.2 


64.8 


80.0 


over. 


10,000 


10,000 


10,000 


10,000 


10,000 












stag. 


10,000 


10,000 


10,000 


10,000 


10,000 


10,000 


10,000 


10,000 


10,000 


10,000 



Table 2: Survey of matched T = L/4 lattices, with statistics for overlap and staggered fermions. 
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Figure 1: Low-energy spectrum of the unfiltered and two filtered versions of the Dirac operators 
D st and D ov [after chiral rotation as in pHjl ] on three typical configurations at /3 — 7.2 and on 
a selected one where the overlap charge depends on the smearing level (rightmost panel). 



eigenvalues form near-degenerate pairs which sit close to an individual overlap eigenvalue. In 
particular the right number of "would-be" zero modes separates from the rest of the spectrum 
and clusters near the origin, namely 0,2,8 for q — 0, 1,4, respectively. This is a manifestation 
of the (approximate) index theorem for staggered fermions which holds both in 2D [7 j and 
in 4D jHl El EI] ■ However, the similarity extends to the higher modes, and one can take the 
"fingerprint" of a configuration likewise with D ov or D st , the two differ just by a trivial rescaling 
factor and the two-fold degeneracy in the latter case. This is the basis of the rooting procedure 
for staggered fermions, and it is thus interesting to perform, in a first step, a scaling analysis 
for some observables which can be formed from the eigenvalues only. 



3 Scalar condensate 

For overlap fermions, the (bare) scalar condensate is unambiguously defined as [21] 

Denoting the eigenvalues of the massless overlap Dirac operator by A, and remembering that 
we work with p — 1 , the reweighted condensate is 

y ./a (det(D° v ) Jv fE n ^(L ) m 

Asca VH N v ml ^ (1— m/2)\+m' j , / n ov\ TT/Yi \\ i \ /1 o\ 

— = U (det(D%m ' det ^) = n((l-y)A+-) (12) 
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where the sum runs over the full spectrum. These eigenvalues occur either in complex conjugate 
pairs or as isolated chiral (doubler) modes at A = (2). Finally, one can rewrite (|T2*j) as 



e L 2 (det( J D^) JV /) 



A = (A- 1 -2- 1 ) 1 (13) 



where A is purely imaginary and the primed sum excludes the doubler modes at A = 2. 

In the staggered case we follow [22] and implement the (bare) 1-flavor condensate through 

*f = ~YT^ xx) ' (14) 

where the purpose of the factor 1/2 is to compensate the two-fold degeneracy of the staggered 
formulation in 2D. Denoting the eigenvalues of the massless staggered Dirac operator by A (they 
show up in complex conjugate pairs with zero real part), the reweighted condensate is 



e 2L 2 (det(D%) N f/ 2 ) 



detK) = n( A+m ) • (15) 



with the sum and product running over the entire spectrum. 

Finally, there is no renormalization factor to be taken into account in a 2D theory with cut- 
off effects of order 0(a 2 ), if one insists on a massless renormalization scheme. This follows from 
the general expansion Z = l-+0(e 2 ) which, due to the dimensionful coupling, reads Z = l-+0(a 2 e 2 ) 
and thus tells us that there is no way to separate a Z-factor from intrinsic 0(a 2 ) effects. 

An example of the condensates ()13(ll5j) with Nf = 1 in one of our zero temperature geometries 
is shown in Fig. |21 The overlap yields a smooth curve which, for m— >0, is consistent with the 
Schwinger value (JHJ), marked by an asterisk. This holds both with a standard Wilson kernel 
and with the filtered version. By contrast, the staggered condensate (at any filtering level and 
any (3) tends to zero, if the chiral limit is performed at fixed lattice spacing, and taking, in a 
second step, the continuum limit will not change this, thus 

Jim Jim X ^aKa) = Q (16) 

a— >0 e 

Given the dramatic difference between the two staggered curves in Fig. |2J one wonders whether 
the staggered answer could be useful, if one uses only the data above some m min to extrapolate 
to the chiral limit. Obviously, there is no canonical definition of such a m m ; n , but the staggered 
results in Fig. El suggest that it might be less ambiguous for the filtered operator. The question 
can be formalized by asking whether it is possible to reproduce the Schwinger result (JHJ) with 
rooted staggered fermions, if one considers the reverse order of limits, i.e. whether 

m-tOa-tO e 27T J / 2 

The goal is to show that (only) the second order of limits works with staggered fermions, while 
the overlap discretization is correct with either order. 

It turns out that eqn. (JTTJ) together with our definitions for Xsc/ St does not make sense, 
neither for staggered nor for overlap fermions. This is because the definitions (jllj) and (JHJ), 
when evaluated with a positive quark-mass, lead to a logarithmic divergence in the cut-off. 
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overlap, no smearing 
overlap, 1 APE step 
staggered, no smearing 
staggered, 1 APE step 



0.15 
m/e 



0.3 



Figure 2: Bare overlap and rooted staggered Nf — 1 condensate Xsca/e at 0=7.2, plotted versus 
the quark mass. Here and in subsequent figures, three lines indicate a lex-error band. The 
overlap result changes little, if a UV filtered Wilson operator is used instead of an unfiltered 
kernel. By contrast, in the staggered case this makes a big difference - only the filtered variety 
shows a clear separation into a regime (here: m/e > 0.02) where the staggered answer is 
meaningful, and a regime (m/e < 0.02) where lattice artefacts overwhelm. 



The logarithmic divergence shows up only in the condensate at finite m, with a coefficient 
that vanishes as the quark mass tends to zero. That such a "soft breaking" through mlog(A 2 ) 
terms must happen is obvious from the free case, where elementary manipulations yield 



m J p 2 + m 2 Jo Jo A 2 + m 2 
A „ Am , /A 2 + m 2 - 



27r lim / dX — = lim 7rmlog(- 

A^oo Jo X 2 + 111 2 A^oo V 



m 2 



Fig. El collects our Nf = l overlap data at zero and finite temperature, evaluated at m = 0. 
There is no sign of a logarithmic divergence, and we get an acceptable fit if we use 

xUm=0,l/M =A + Bj(3 (ig) 

with a common A and slope parameters B n which account for the 0(a 2 ) effects at filtering 
level n. At zero temperature our continuum result A = 0.1629(47) is compatible with the 
Schwinger value (JHJ). The result A = 0.0700(16) from our thermal geometries indicates that at 
this temperature the massless continuum condensate is considerably reduced. The main lesson 
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Figure 3: Bare overlap condensate X°ca/ e with m = at zero temperature (left) or finite 
temperature (right) versus (ae) 2 . The three filtering levels have a universal continuum limit, 
and the extrapolation works well. The zero temperature result is consistent with (jHJ). 
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Figure 4: Bare Nf = l condensate Xsca/e with overlap (left) or rooted staggered (right) fermions 
at zero temperature and m/e = 0.1 versus (ae) 2 . The three filtering levels have a common 
logarithmic divergence. Adopting the ansatz ([20)1 . the divergence may be subtracted for either 
discretization, but the so-defined continuum limit is non- universal. 

is that with overlap fermions one can indeed take the chiral limit first, letting a^Oina second 
step, and gets the correct answer. 

Fig. |U shows the Nf — 1 zero temperature overlap and staggered condensates, evaluated at 
m/e = 0.1. Now, the divergence in the cut-off is clearly visible. Indeed, we get an acceptable 
result from a correlated fit to all data (in one formulation) at fixed m/e with the ansatz 

X Um/e,l/(3,n) =A+ +C]og (i/0) + Mf^+^V^ + WtfV/* 1 (20) 
e 

where all coefficients depend on m/e, and A, C are common to all filtering levels, while B n , 
£)(3/2,2) depend on the level n. Thus, besides the leading 0(a 2 ) cut-off effects for each filtering 
level, we include 0(a 3 ) terms in levels 1,0, and 0(a 4 ) terms with no smearing. With rooted 
staggered quarks everything is analogous, i.e. the ansatz (}20| works again. The fine dotted lines 
beneath the data indicate what remains if one subtracts the so-determined logs plus higher order 
terms and stays with the constant plus linear part. For each m/e > there is a well-defined 
continuum limit with such a procedure, but there is no reason to expect it to be universal. 



The difference between the "minimally subtracted" overlap and staggered condensates that 
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Figure 5: Left: Free condensates versus 1/(3 (in our zero temperature geometries, extended 
towards larger /3), at fixed m/e = 0.1. The fit to eqn. (j2IJ) uses the data below 0.4, i.e. L>16. 
Right: Difference of the free staggered and overlap condensate versus 1/(3. The data to the left 
of 0.4 are fitted with a 0(1/ (3) + 0(1/ (3 3 / 2 ) ansatz, shown together with its leading part. 



survives in the continuum can be traced back to the difference between the two free condensates. 
The latter difference has a finite continuum limit, since each of the free condensates contains the 
same logarithmic divergence. This point is illustrated in Fig. El On the left the free condensates 
X°re{ st /e at fixed physical quark mass are plotted versus 1/(3 (here, we refer to Tab. [TJ i.e. (3 
merely encodes for the geometry). The curves have a common log term; they fit to the form 



ov/st 
Xfree 



m/e,l/(3) 



A°:J e st (m/e) + B°;J c st (m/e)/P + C frce (m/e) log(l//3) 



(21) 



where A {Tee , B {ree depend on the discretization, while Cf ree is universal. The same conclusion 
is reached via directly plotting the difference, as shown on the right; the curvature decreases 
towards the continuum, and there is no sign of a remnant log. 

The lesson is that one should define the massive continuum condensate in 2D through 



Xsubt(W e ) 



ii A 2 + m 2 



Pfrcc (A) 



(22) 



in infinite volume and regard its computation with overlap and staggered fermions a direct test 
of universality. Here, the underlying assumption is that asymptotically the eigenvalue density 
p(X) agrees with the free one [from a glimpse at (fTSj) . one learns that Pf ree (A) — > ir\]. 

A legitimate strategy to implement (}2*2*|) on the lattice is to first subtract the free case in 
the same formulation and geometry, followed by an extrapolation to infinite volume. The first 
step may be done analytically, the corresponding expressions read 



£ Tr/^M-^l- 1^(0); 

(p-™) (l+ D ^- p) )+m 
2 + M-£c, + i]r 



(23) 



1 L/2-1 
^ fci,fc 2 =0 



2m 



m 2 + sj + s 2 



(24) 



8 




Figure 6: The "naive" subtracted condensates (J251I2E1) at (3 — 7.2, together with x sca and x& 
for overlap (left) and staggered (right) fermions. Note the power-like IR-divergence in x°ubf • 




m/e m/e 



Figure 7: The "primed" subtracted condensates (jSHllSDI) at (3 = 7.2. No IR-divergence is left. 



with c M = cos(^/c M ) and = sin(^fc Ai ). In the ("naive") subtracted condensates 

xZu(m/e,l/P,V) = &(m/e,l/P,V)-j<g aB (jn/e,V) (25) 
X f uht (m/e,l/(3,V) = X f ca He,l/(3,V)-xt e (m/e,V) (26) 

the limit V — > oo must be taken before m — > 0, since the free case contains a massless mode, 
and this reflects itself in the behavior X°ubtV e— *~ 2/(emV) for a small mass. 

Looking at ()23[ E3| the piece that makes the chiral limit in a finite volume singular is easily 
identified, it is the term with k\ = k2 = 0. Hence, a better choice is to subtract the free case in 
the same geometry without the non-topological zero-modes . With 

xSW = ^£'&{%M- 1 (i-^fl£.(o))j (27) 

= z^'i^rfm < 28 > 

[the primed sum skips the (0,0) contribution, otherwise the summation is as in (|23 | 124 )) ]. one 
may define the "primed" subtracted condensates 

X ' s :u(m/e,l/(3,V) = X Z> (m/e, 1/(3, V) - xg(m/e, V) (29) 
X^ bt (m/e,l//3,V) = x s L{mle,ll(3,V)-x!^{m/e,V) (30) 
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m/e m/e 



Figure 8: The "primed" subtracted condensates (j291l3U|l for overlap (left) and staggered (right) 
fermions with 1 smearing step. Note the small scaling violations at any m/e with overlap 
fermions, while they get progressively worse for small m/e in the staggered case. 



that avoid the IR singularity in (f23| l2T)|) . In other words, (J2*H| l3T)|) is the only lattice imple- 
mentation of (j22| which is simultaneously UV and IR finite and thus permits to skip the limit 
V ^>oo, because finite volume corrections are (asymptotically) exponentially small. Note that 
this procedure differs from a standard renormalization in QCD. Since the theory is super- 
renormalizable, there is no bare parameter to be adjusted and no counterterm to be added to 
the Lagrangian, the only option is a pure vacuum reordering. 

The "naive" and "primed" condensates (|25[ l27)|) and (|29| l3T)|) are shown for one coupling in 
Fig. Hand Fig.[7[ both for overlap and staggered quarks. The former version clearly exhibits the 
IR divergence X°^it( m / e ~ >®)l e ~ * — 2/(emV), due to the subtraction term in (|2*51 EE} - In the 
latter version, the leading finite- volume effects are tied to the lightest particle in the interacting 
theory, the Schwinger particle with mass ej '^/TT+0(m), and thus exponentially small. 

Having defined the "primed" subtracted condensates, we are in a position to let, for any 
fixed m/e>0, the lattice spacing tend to zero 

lim oV ^ X ^ t (m/e,l//?) (31) 

lim^;iKe,l/^ (32) 

and Fig. |H1 indicates that there is, at this stage, a practical difference among the two formula- 
tions. In the overlap case cut-off effects are mild, and this means that the continuum limit can 
be taken, from the /3-values available, over the full mass range shown. With rooted staggered 
quarks, scaling violations get progressively worse for small m/e. 

Fig. El indicates that, once the fermion mass is large enough to be in the staggered scaling 
regime, both formulations agree in the continuum, and the agreement seems not to be tied to 
a particular Nf. Comparing our [Nf = 1) continuum result Xcont( m / e = 0-1) = 0.142(4) to the 
values 0.1171 by Hosotani [2E] and 0.1277 by Adam |2Z| we find a disagreement at the 6.2a 
or 3.6a level, respectively. This might indicate a finite-volume effect or a limited precision of 
the latter calculations. Since our boxlength is ~ 5 times larger than the correlation length of 
the lightest particle in the massless Nf = 1 theory, and in this regime finite-volume effects are 
exponentially small we consider the first option with our data unlikely. 

As a last step, one my now take the chiral limit. In the staggered case we find that after 
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Figure 9: Continuum limit, via (f2T)J) without the log (full lines), of Xsubt/ e a ^ fixed m/e — 0.1 with 
overlap (left) and staggered (right) quarks for Nf — 0, 1,2 (top, center, bottom). Comparing to 
the strictly linear extrapolation (dotted lines) gives a measure for the systematic error. 



the a— >0 limit has been taken with the ansatz 

X ' suh t(m/e,l/P,n) =A+ + Q + /2 

e 

we have a continuum curve that covers the range 0.05 < m/e < 0.15, outside the ansatz ()33|) 
yields an unacceptable y 2 . Using this segment and a linear ansatz for the extrapolation m^0 
we find that the massless staggered Nf — 1 condensate agrees with the analytical result (jHJ). In 
the overlap case the continuum curve would extend to m/e = 0, but as a consistency check we 
perform the same extrapolation, getting again a result compatible with (jHJ). Both extrapolations 
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Figure 10: Chiral extrapolation of the continuum extrapolated primed condensates (jHH E2|)- 
In the staggered case the continuum limit can only be taken for m/e > 0.05 from our /3-values, 
and with overlap quarks the same cut is applied for consistency. A linear extrapolation m — >0 
agrees with the analytical result (jHJ) for either discretization. The error-band is only statistical. 



are shown in Fig. Thus, for overlap fermions our data support the universal behavior 



lim lim 

a— >0 m— >0 



Xsubt(We,a 2 



lim lim 

m— >0 a— >0 



Xsubt(We,a 2 ) exp( 7 ) 



e m— »U a^U e 27T 3 / 2 

For staggered fermions, one the other hand, the suggested non-commutativity phenomenon 



(34) 



lim lim 

a^0 m— >0 



Xsubt(™/e,a 2 







lim lim 

m— >0 a— >0 



Xsut,t(We,a) exp( 7 ) 



e nnUa-»U e 27T 3 / 2 

means that rooted staggered fermions can be used to reproduce the Schwinger result 
only if the right order of limits is chosen. 



(35) 
), but 



4 Topological susceptibility 

Another interesting observable to study the effects of dynamical fermions is the topological 
susceptibility which, in the context of this note, shall be defined (in the continuum) through 

(det(D+m) N f q 2 ) 
X ^ = ^V(det(D +m rf) ■ (36) 

The main difference to the scalar condensate is that the topological susceptibility depends only 
on the sea-quarks, thus offering a potentially cleaner view at the effects of square-rooting the 
staggered determinant to get Nf = l. 
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For staggered quarks, the definition (J36|) . taken in fixed volume, reduces to 

Xgp g (det(^)^// 2 g 2 ) 

e 2 L 2 (det( J D^) JV // 2 ) ' 1 ' 

while for overlap quarks, the implementation is 

xSp g (dfit(£%)"/ g 2 ) 

e 2 L 2 (det(D£)*/) ' 1 J 

where det(-D^) and det(D^) are defined in ([TBI [T2J. The sum is over the full spectrum of the 
massless operator, and we apply the overlap definition (JTUJl of the charge q in either case. 

Let us begin with a check that a reasonable number of our zero temperature (3- values is in 
the scaling regime before reweighting. Fig. [TT] contains our quenched topological susceptibility 
data, and /3 > 3.2 seems sufficient to be in the regime with only 0(a 2 ) effects. 

The first test is whether overlap and rooted staggered quarks yield the same result in the 
continuum, if the topological susceptibility Xtop _1 ( m / e )/ e2 is taken at a fixed non- vanishing 
quark mass. Fig. El presents the outcome for m/e = 0.1. In either case the three smearing 
levels seem to have a common continuum limit and this is why we adopt a correlated linear 
fit. The staggered continuum value 0.00958(42) is consistent with 0.01011(18) from all overlap 
data. This is of course not a proof but good numerical evidence that rooted staggered quarks 
yield the correct continuum limit for the topological susceptibility at finite quark mass. 



0.038 




1/P 

Figure 11: The zero temperature quenched topological susceptibility Xtop / e<2 versus (ae) 2 . 
The symbols for nsmear = 0,3 are slightly offset for better visibility. Note the steep slope - 
in spite of being in the scaling regime for (3 > 3.2, cut-off effects can be large. At /3 — 9.8 the 
charge q is independent of the smearing level for all our configurations. 
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Figure 12: At the fixed quark mass rn/e — 0.1 a correlated fit to all staggered Nf — 1 data (top) 
suggest a continuum limit 0.00958(42) which is consistent with 0.01011(17) from all overlap 
data (bottom). In the staggered case, the first smearing renders the extrapolation much flatter, 
while the overlap has a flat behavior with any level of filtering (note the difference in scale). 
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Figure 13: The remnant Xto P / e2 a ^ Ta/e — Q versus in log-log scale. This staggered artefact 
seems to disappear with a power-law in a, for Nf = l (left) and Nf = 2 (right), but the power is 2 
only for smearing level 0, while it is 4... 8 for higher filtering levels (over our range of j3- values). 

The second test concerns the chiral limit, where x°o P (0) = ( as i n the continuum [28]). while 
Xto P (0)>0 is a pure discretization effect. Fig. HH1 shows the remnant staggered susceptibility at 
m/e = versus (ae) 2 , in a log-log representation. With unfiltered staggered sea-quarks it seems 
to disappear in proportion to a 2 , regardless of Nf. For the filtered variety we find a slope ~2 
(i.e. dominating 0(a 4 ) cut-off effects) over the range of accessible couplings for Nf = l and even 
~4 (i.e. dominating 0(a 8 ) effects) for Nj = 2. Still, it is conceivable that this slope eventually 
flattens out and the asymptotic cut-off effects might be 0(a 2 ) with any level of filtering. This 
would mean that smearing renders the coefficient in front of the 0(a 2 ) term so small that the 
"subleading" 0(a A ) or 0(a 8 ) terms would numerically dominate over a substantial range of 
couplings. In any case, the important news is that the non-commutativity phenomenon in the 
chiral condensate is not replicated here - for the topological susceptibility the staggered answer 
is correct even if the chiral limit is taken before the continuum limit. 



5 Partition function and Leutwyler-Smilga sum rules 

In the e-regime of QCD [2E] (m£V<Cl, but still with a large box-size, i.e. L^>l/(2F 7r )~lfm, 
c.f. the discussion in |29j ) the log of the partition function is known analytically [28J 



Z ( m ) = e mEycos W (39) 

with sub-leading corrections of order {mLV) 2 . The partition function in a sector of fixed 
topological charge q follows by taking the Fourier transform 



Z q {m) = J de e iqd Z g (m) = I\ q \(mY,V) (40) 

where I\ q \ is the modified Bessel function. 

By differentiating with respect to the quark mass and setting the latter to zero, Leutwyler 
and Smilga obtained the first sum rule for the inverse eigenvalues of the massless Dirac operator 

(^» (imag(A n )Sy)2) g = 2(\q\+N f ) (41) 

where the primes indicate that in a sector with topological charge q ^ the sum is over the full 
spectrum without the |g| zero modes and the importance sampling is to be carried out with the 
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Figure 14: The Nf — 1 partition function Z\ q \ (normalized by the sum over all q) at (3 = 3.2 
without filtering (top) and at (3 = 9.8 with 1 smearing step (bottom) versus x = mY>V. For 
sufficiently small mass the overlap version (filled symbols) is always a decreasing function of 
\q\, as in the continuum. The rooted staggered version (crosses) joins, for large mass, the 
associated overlap Z\ q \/Ze =0 . The full lines represent the parameter-free prediction l40 |) . 



16 



2 
1.8 
1.6 
1.4 

1.2 

s 

W c 1 



P 



0.8 
^ 0.6 

0.4 
0.2 


2 
1.8 
1.6 
^ 1.4 

> 
W 

c 

£ 

® 0.8 
0.6 
0.4 
0.2 








0.2 



1/P 



- 




i 




i 


1 




i 


- 


- 










• 

X 


q=0, overlap 
q=0, staggered 


- 


- 












■ 

+ 


lql=l, overlap 
lql=l, staggered 


- 


— 








X 






lql=2, overlap 


— 












* 


lql=2, staggered 




— 
— 






X 










— 
— 


— 
















— 








• 




+ 








_ 















* 


_ 


— 




— *- 


H 


-* 








— 






1 




1 


1 




i 













0.2 


0.4 












1 




1 


1 




i 






• 


q=0, overlap 














X 


q=0, staj 


jgered 












- 
- 


■ 

+ 


lql=l, overlap 
[ql=l, staggered 










- 
- 






lql=2, overlap 
lql=2, staggered 








X 




— 

— 
















— 

— 






• 


* 








• 




_ 




t- 




i 

T 


+ 




+ 


_ 


— 








-■- 






* 


— 






i 




1 


1 




i 





0.4 



Figure 15: The Leutwyler-Smilga sum rule (|4"2j) for Nf = 1 active quarks versus (ae) 2 , at 
smearing level (top) and 1 (bottom). Note the dramatic difference in the scaling behavior 
between overlap and rooted staggered fermions. The deviation, in the continuum, from the 
Leutwyler-Smilga prediction (full lines) is likely a finite-volume effect. 
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"primed determinant", i.e. without the zero modes. In 4D the l.h.s. is quadratically divergent, 
since the eigenmode distribution is asymptotically proportional to A 3 . In 2D the divergence is 
logarithmic, since the eigenmode distribution is asymptotically linear in A [cf. (I22j) ]. One way 
out is to subtract the free case in the same volume. An other one is to consider the difference 
between two topological sectors, since the UV-behavior is independent of the charge. 

We start with the conjecture that eqn. (}39|) and thus (J4T^I4T| hold in the 1-flavor Schwinger 
model, too, albeit with re- interpreting E as the analytically known 1-flavor condensate (JHJ). 

In Fig. 1141 we show the Nf = l partition function on a coarse and a fine lattice, without and 
with one APE step, respectively. In the first case, one sees a significant difference between the 
overlap and the rooted staggered answer for small quark masses. For larger j3 and with just one 
smearing step this difference disappears and we did not find any indication for a disagreement 
in the continuum limit at any m/e. For x<Cl the continuum extrapolated data agree well with 
(2 — 5 q fl)I\ q \(x)/ exp(x) which, with our choice for E, is a parameter-free prediction. We like to 
point out that this figure is reminiscent of the situation in QCD, see Fig. 3 of Ref. [3*0] . 

To verify ()41|) we subtract the free case, i.e. we check whether 

(imag(A n ) EV) 2 ~ ^« (imag(A^ e ) ZV) 2 \ = 2(\q\+N f ) ^ 

with A ra denoting the (purely imaginary) n-th staggered or chirally rotated overlap eigenvalue 
(i.e. the A of (HHJ)), and in the staggered case an additional factor 1/2 on the l.h.s. is needed (or 
1/4 in 4D). The prime on the sum indicates that on topologically nontrivial configurations the 
|g| true overlap zero modes or the \q\ "would be" zero modes on either side of the staggered 
spectrum are excluded. Since also the free overlap and staggered Dirac operators have 2 and 4 
(non-topological) zero modes, respectively, the A^ ee denote the free eigenvalues without these 
zeros. Thus the sum over n does not include the same number of terms in the free and interacting 
case. The prime above the expectation value indicates that in a q^O sector the \q\ zeros are 
removed from the determinant. Fig. contains our results for the sum rule (|4*2*|) . where our 
choice for the interpretation of E leads again to a parameter-free prediction. The deviation 
from the Leutwyler-Smilga value is likely a finite-volume effect. In any case, there is no sign of 
a difference between the rooted staggered and the overlap answer in the continuum. 



6 Heavy quark potential 

The heavy-quark (HQ) potential is an interesting observable, since it is easy to measure in the 
Schwinger model and the effect of dynamical fermions is clearly visible. We define the HQ 
potential V(r) in a finite volume L x T via the Polyakov loop correlator 

e -V(r)T = (p( )pt( r )) . (43) 

In the quenched case, and with T — > oo, it has a simple saw tooth shape, i.e. it starts at 
zero, raises linearly for < r < L/2, and then it decreases linearly, until it reaches at r = L 
again. This follows from the fact that the quenched Wilson loop in infinite volume satisfies an 
exact area law WN f=0 (r,t) = (i 1 (/3)/i (/3)) rt [in lattice units], thus the force in physical units is 

F Nf=0 ,T=oo(r) (h{(3) 



e 2 



"^W- (44) 



In the Nf = 1 case an exact solution on the torus is known for m = J jT| • but - as far as 
we know - not for massive fermions. However, we expect the screening behavior found in the 
massless case to also apply, on a qualitative level, for a small enough fermion mass. 
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Figure 16: The quenched zero temperature HQ-potential is a saw tooth. After reweighting with 
staggered (top, rooted for Nf — 1) or overlap (bottom) fermions (m/e — 0.1), screening is seen. 
Throughout, the same smearing is applied in the HQ line and in the fermion operator. Except 
near the endpoints, the improved HQ action leads to a mere downwards shift of the potential. 
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Figure 17: Top: With our values, the quenched heavy-quark force F/e 2 at r = L/4 is in 
the scaling regime. For large enough /3, the three HQ actions agree in their central values and 
errors. Bottom: The Nf = l HQ force F j e 2 at r = L/4 with overlap or rooted staggered fermions 
(filtering level 1 applied in both and the HQ action) seems to have a universal continuum limit. 
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Figure 18: The difference between the HQ potentials with smearing-level and 1 (at r = L/4) 
versus ae, in log-log scale. It disappears in proportion to a - the offset line has slope 1. 



In Fig. we show an example of our HQ potentials {(3 = 12.8, lattice units). In this 
regime the quenched prediction ()44j) works very well, except that the finite temperature causes 
a flattening of the peak region. If the Polyakov loops are built from smeared links instead of the 
original ones ("improved HQ action", see [HHIH2]), the potential is merely shifted downwards, 
except near r = 0, L where some curvature is introduced. In other words, the force is unaffected 
by such a modification of the HQ action, provided it is measured at a distance larger than the 
number of smearing steps applied. Furthermore, the reduced HQ selfenergy is supposed to lead 
to a smaller error bar (originally [31] . see jH2] for details), but here it seems unaffected - we will 
come back to this point. After reweighting to Nf — 1 or Nf — 2, screening is seen (as expected), 
yet there is no visible difference between staggered (top) or overlap (bottom) dynamical flavors. 
Here we show only the "diagonal" data, where the same smearing level is applied in the HQ 
action as in the fermion operators, but the "offdiagonal" data are similar. 

For a detailed analysis we shall concentrate on the force at a given (large enough) distance, 
and we choose to present a scaling analysis for F Nf=1 (r = L/4)/e 2 . Before doing so, it is worth- 
while to check that our /3-values are such that the quenched counterpart, F Nf=0 (r = L/4)/e 2 , 
is in the scaling regime, and the top of Fig. IT7I shows that this is indeed the case - regardless 
of the HQ action (smearing level of the Polyakov loop links) used. When employing overlap 
or rooted staggered fermions to have one active flavor with m/e = 0.1, nothing fundamental 
changes, as shown in the bottom of Fig. El in fact, at the 4 couplings where we have data 
with either discretization, the results are fully compatible already at finite (3. Thus, in the case 
of the HQ force a rooted staggered field seems to yield the same continuum limit as a single 
overlap flavor. 

Let us finally discuss an interesting aside. It has been argued (originally see [32| for 
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details) that the selfenergy fa of a quark in the static approximation is directly related to the 
noise in an observable with such a quark. The important point is that in 4D the selfenergy in 
physical units diverges linearly in the inverse lattice spacing [3*H] 

8m = const /a [4D] . (45) 

If changing the discretization amounts to a replacement const' — > const", then smearing becomes 
more important on fine lattices. What we wish to point out is that the situation in 2D is just 
opposite. Here, the selfenergy vanishes in proportion to the lattice spacing 

6m = A 2 m a [2D]. (46) 

As a consequence, the usefulness of smearing decreases in 2D towards the continuum. This is 
clearly seen in the top of Fig. ^| In the rightmost point (/? = 3.2) the smeared HQ action leads 
to a smaller error bar, while this effect quickly disappears towards the left. This is matched by 
the behavior of the shift of the HQ potential [in physical units], brought by a single smearing 
step, as shown in Fig. HH1 Such a pattern is expected, if the first smearing step amounts to a 
replacement A H q^A H q, and it is reassuring to see that this rule holds for any Nf. 



7 Determinant ratios 

The last topic that we wish to discuss is how well the rooted staggered determinant manages 
to approximate the (one flavor) overlap determinant. It has been shown - both in 2D |7] and 
in 4D jH] - that the low energy eigenvalues eventually coincide (apart from an overall rescaling 
factor and the 2 d//2 -fold degeneracy), in the limit of fine lattice spacings. However, in the UV 
region the two spectra are very different, and, if the resolution is made finer in a fixed physical 
volume, the total number of modes grows in proportion to l/a d . Therefore, it is not clear 
whether the full one-flavor determinants would eventually coincide. In fact, det(.D ov ) and the 
2<V2_foid root of det(£> st ) are never equal (typically, they differ by many orders of magnitude), 
but the right question to ask is whether the two formulations give the same answer - modulo 
cut-off effects - for the determinant ratio on two arbitrary configurations, i.e. whether 



det(D m (U)) 



2 d/2 



det(D st (U)) 



0{a 2 )) . (47) 



det(D ov (U')) " \ det(D st (U')) 

We are in an excellent position to test (|47j). since we have all eigenvalues with either dis- 
cretization in a number of lattices with fixed physical volume (see Tab. [Q). Introducing 



b ^ U > U > ~ 2^ [ ° g {det(D%(U>))) 



(49) 



the goal is to show that the unfiltered ratios differ by cut-off effects only, relation (}4Tj) or 

S%(U,U') = S*(U,U') + 0(a l ). (50) 
For this it is sufficient to show that an analogous relation holds among the filtered operators 

[/'«) = S%{UM, t/' (1) ) + 0(a 2 ) , (51) 
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Figure 19: Log-log plot of the rooted staggered versus the overlap determinant, at fixed physical 
mass m/e — 0.1 and /3 = 1.8, 5.0, 9.8 (from top to bottom), without (leftmost column) and after 
one step of filtering (second column). In addition, we plot smearing level 1 versus level for 
the overlap (third) and the rooted staggered determinant (fourth column), which are known to 
deviate through 0(a 2 ) effects. Throughout, the black line is the identity, not a fit. 



where the superscript (1) refers to the smearing level 1, since we already know that the left-hand 
sides and the right-hand sides of (joTH 15*1]) satisfy 

S%(U,U') = S%(UW,UM) + 0{a 2 ) (52) 
S%(U,U') = S^(U (1 \U^) + 0(a 2 ) (53) 

respectively, because filtering amounts to an 0(a 2 ) redefinition of the operator. 

The leftmost column of Fig. contains scatter plots of the two sides of ()50|) in a fixed 
physical volume and at a fixed quark mass m/e. Each dot represents a configuration U, while 
U' is an artificial background which realizes the ensemble mean (over U) of (|4"%1 l4T)|) . Obviously, 
the correlation on an arbitrary background U gets tighter and the pertinent slope moves closer 
to 1 as the lattice spacing decreases. Note that this holds regardless of the topological charge, 
since our ensembles contain a Gaussian distribution of q at each {3. The second column shows 
the scatter plots of the two sides of (}5T| . where the correlation is much better than before. 
Finally, the third and fourth columns contain the scatter plots for the known relations (|52|) 
and (|53jl . Qualitatively, they look very similar to the first two columns, and this gives some 
confidence that (JSTJl l4"Tj) might actually hold true. In any case, the interesting news is that the 
1-filtered staggered action generates, on a fine enough lattice, a distribution that is closer to 
the the 1-filtered overlap one, than the latter is to the unfiltered overlap distribution. 

To make the discussion a bit more precise, we introduce a quantity designed to measure the 
deviation of the correlation from the identity. We use e.g. for the first column 

X 2 = ^ 7 yE(^og(det( J D st ))-log'(det( J D ov ))) 2 (54) 
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Figure 20: Scaling of x 2 > defined in (|55|) . versus (ae) 2 . For the diamonds/bullets it is known 
that they represent an 0(a 2 ) redefinition of the operator, for the crosses/stars it is conjectured, 
because of the striking similarity of the fall-off pattern. 



where log'(det) = log(det) - l°g(det) 

W{ u ) u> 

and expect that it decreases, in a fixed physical volume, as the lattice spacing gets smaller. Of 
course, this y 2 is not a physical observable, and one should not expect it to vanish in proportion 
to 1/ ' (3 = a 2 . Still, the data in Fig. |2"U1 suggest that it falls of - roughly - with a power law in a. 
We tested other m/e which are in the scaling regime with our /3-values and found qualitatively 
the same behavior. The important point is that there is no obvious difference among the 
pattern of the four coefficients pertinent to the four columns of Fig. El Only the prefactor is 
different, telling us that the correlation among the two formulations at smearing level 1 is better 
than the internal overlap correlation between level and 1, and the latter is better than the 
internal staggered and the staggered-to-overlap correlation without smearing. In other words, 
the scaling of the coefficient (Jo7)|) between the overlap and the rooted staggered contribution to 
the effective action (normalized per physical volume) behaves in the same way as the one for 
the internal overlap or internal rooted staggered correlation between two smearing levels. Since 
the latter are known to deviate through 0(a 2 ) effects, we rate this as strong evidence that also 
the overlap and rooted staggered determinants differ by 0(a 2 ) terms, cf. (|3Ul I5T]) or (jUJ). 



8 Summary and discussion 

We have attempted scaling tests for dynamical overlap and rooted staggered fermions. The 
Schwinger model was chosen, because it permits to reach high numerical accuracy while the 
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conceptual issue of square- rooting the staggered determinant remains for Nf = l. 

First, we have looked at the scalar condensate at fixed physical quark mass. For m/e>0 the 
bare staggered and overlap condensates (jlH IT4*|) diverge logarithmically in 1/a. Subtracting the 
free case without its non-topological zero modes defines a UV and IR safe observable at finite 
quark mass, and Fig. El suggests that it has the same continuum value for staggered and overlap 
fermions. Hence the analytic value (jHJ) for the condensate in the massless 1-flavor theory can 
be reproduced with rooted staggered fermions, if the limit m— >0 is taken after letting a—>0. 
In the overlap formulation the two limits may be interchanged, while in the staggered case the 
reverse order of limits yields an exact (but incorrect) zero. We consider the non-commutativity 
phenomenon (jHo^) important, because the Schwinger value (JHj) reflects the axial anomaly, and 
the second version indicates that rooted staggered fermions do get the anomaly right, if the 
continuum limit is taken at finite quark mass. 

The second observable has been the topological susceptibility. Both in the continuum and 
with overlap sea-quarks it vanishes at zero quark mass. Hence Xtop( m = 0) is a pure lattice 
artefact, and it seems to disappear like a 2 for an unfiltered ("naive") staggered Dirac operator, 
while in a filtered version 0(a 4 ...a 8 ) terms would numerically dominate over a wide range of 
couplings. On the other hand, with a finite quark mass, there is no obvious difference among 
the two formulations, and we have verified that one sees pure 0(a 2 ) scaling, as expected. Fur- 
thermore, the continuum limit of the topological susceptibility with overlap or rooted staggered 
fermions at fixed physical quark mass is consistent within errors. 

Our third observable, the partition function, is a bit more refined, but similar in spirit, to 
the topological susceptibility. Again we find, at fixed m/e, a continuum value, for arbitrary 
charge q, that is consistent within errors. Also the fourth test, a check of the Leutwyler-Smilga 
sum rule ()42j) has not revealed any difference of the two formulations in the continuum limit. 

The finite temperature heavy-quark force at fixed physical distance has been our fifth test. 
It is interesting that in the Schwinger model the effects of virtual quark loops are so pronounced 
and so easy to measure. We find no deviation between the two fermion discretizations, since at 
filtering level 1 overlap and rooted staggered sea-quarks agree already at the /^-values considered. 
It is noteworthy that the a-dependence of the heavy-quark selfenergy in 2D is entirely different 
from the one in 4D. As a consequence, a fuzzed heavy-quark action increases the signal-to- noise 
ratio only on rough lattices in 2D, while this trick helps only on fine lattices in 4D. 

Finally, we have investigated the logarithmic determinant ratio of overlap and rooted stag- 
gered fermions as a function of the lattice spacing in a fixed physical volume. We find strong 
evidence that, at a given filtering level, this ratio is constant over the entire configuration space 
up to 0(a 2 ) effects, since the pertinent \ 2 disappears in the same way as for (j^SlESl) which are 
established relations. Of course, this is a numerical argument and not a proof. 

In general, we find cut-off effects with dynamical overlap fermions to be smaller than with 
staggered fermions. Close to the chiral limit, the difference may be dramatic. 

As mentioned in the introduction, the conceptual power of a scaling study is one-sided. We 
might have found an observable where the staggered answer turns out wrong in the continuum, 
and this would have been sufficient to prove that the rooted staggered action does not represent 
a legitimate fermion discretization. We have not found such an observable, and the fact that 
even the condensate in the 1-flavor theory is reproduced correctly (under the right order of 
limits) lets one feel more sceptical whether such an observable exists. Still, we emphasize that 
most of our observables involve only sea quarks, and even the finite-mass condensate does not 
probe the complicated flavor/taste structure that staggered fermions have in the valence sector. 
Hence, a scaling study focusing on 7r and rj' properties might be worthwhile. 
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For current dynamical simulations with rooted staggered quarks the implication is two-fold. 
The good news is that, for any finite sea-quark mass, we find a universal continuum limit in a 
dynamical theory. This is just numerical evidence and not a proof, but if overlap and staggered 
quarks would not yield the same result in the continuum, there is no reason why this difference 
would be particularly small, but non-zero, in the massive Schwinger model. The potentially 
worrisome observation is the non-commutativity phenomenon (|35j) . This is critical, because in 
4D physical results are extracted by fitting the data against predictions from staggered chiral 
perturbation theory (SXPT) i.e. the limits a— >0 and m— >0 are taken simultaneously. It 
is not clear to us whether SXPT would accommodate such a non-commutativity (in the full 1 
theory), but maybe standard physical observables in 4D QCD with Nf>2 are not afflicted with 
this problem anyways. 

Finally, we feel that in future investigations the focus should be on the locality issue men- 
tioned in the introduction. Bunk et al. have shown that the fourth root (in 4D) of D st is not 
a local operator [3], but the question is, of course, whether there is any local operator which, 
when raised to the fourth power, reproduces D st . Obviously, this problem needs to be solved 
to lend a solid theoretical basis to a mixed action approach as it has been explored in Ref . [33] . 
Recently, Neuberger has discussed how the operation of taking the fourth root may be cast 
into a local relativistic framework in 6 dimensions |36| . The question that should be asked, in 
our opinion, is whether it is really necessary to have a local operator D that satisfies D A = D st 
exactly, or whether a version with corrections, say 

IT 4 = (D^ + Oia 2 ) (55) 

for the Green's functions, would suffice to guarantee that a "hybrid" formulation with a rooted 
D st in the determinant and D for the valence quarks yields the correct continuum limit 2 . If so, 
our Fig. Uni suggest that such a D might be built via the overlap prescription. 
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